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Abstract 


This report describes the technique developed for shaping the sound spectrum 
of high intensity sound to achieve specialized response characteristics from single 
horns or multiple horn arrays. The design parameter analysis for the technique 
developed includes computations and computer simulation of horn responses using 
the hyperbolic horn theory. This technique provides the ability to minimize any 
undertesting or overtesting of the test article. 


vi 
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Acoustic Spectrum Shaping Utilizing 
Finite Hyperbolic Horn Theory 


I. Introduction 

This study investigates a technique for shaping the 
sound spectrum of high intensity sound with relation to 
better acoustic environmental testing. The computations 
and computer simulations used in this study are based on 
high-intensity sound generator systems. 

The frequency response characteristics of an acoustic 
noise generating system are dependent upon the noise gen- 
erator, acoustic horn (a coupler), and the region for which 
the noise is intended (such as a test chamber, progressive 
wave tube, or free space, etc.). Each of these elements has 
separate frequency response characteristics, therefore, the 
frequency response characteristics of the entire system are 
a function of these characteristics. 

At the present time, developmental contracts through 
other government and NASA centers are developing 
high-intensity sound generators with broader frequency 
response characteristics and greater sound power capa- 
bilities. All these developmental studies deal strictly with 


the sound generator. This particular study is unique in 
that the sound generator remains constant and the sound 
spectrums are varied by specially designed acoustic horns. 

The design for the acoustic horns is based on acoustic 
horn theory for a special family of horns, which are 
referred to as Hyperbolic Horns. These horns have, theo- 
retically, the unique property of a resonant response char- 
acteristic at frequencies just above their cutoff frequency. 

This study includes a complete mathematical analysis 
for hyperbolic horns, and a digital computer program to 
obtain frequency response data for a given single horn 
with particular response characteristics. From this study, 
a technique is developed for designing multiple horn 
arrays. This technique consists of designing each horn of 
the array slightly different from the other horns so that 
the next response function for the array has a smoother 
energy spectrum than any of the individual horns. The 
smoother spectrum is obtained by effectively averaging 
the outputs from all of the individual horns. 
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The techniques developed in this study will allow a 
design parameter analysis for both single horns and mul- 
tiple horn arrays, with optimized response characteristics 
of the single horn (or of an entire multiple horn array). 

The ability to more accurately simulate broadband 
acoustic spectrums for reliable acoustic testing is ex- 
tremely important to minimize the probability of under- 
testing or overtesting the test article. These design 
techniques can also be applied in the audio engineering 
field, i.e., sound reproduction systems could be made 
better by improving the response characteristics of the 
system loudspeakers, both individually and as an array. 

II. Hyperbolic Horn Analysis 

Basic horn theory for a hyperbolic horn of infinite 
length predicts a very sharp frequency cutoff, depending 
on the horn flare constant and on the velocity of sound in 
the medium which for this study is air. For a hyperbolic 
horn of finite length, the response characteristics as a 
function of frequency are dependent upon the following: 

(1) Flare constant N. 

(2) Family parameter T. 

Note: The value T — 1 in hyperbolic horn equa- 
tion theory results in the exponential horn. 

(3) Length L . 

(4) Termination condition z M at the mouth of the horn, 
i.e., whether it opens into free space or an infinite 
progressive wave tube, etc. 

(5) Velocity of sound in the medium (air). 


Expressions for the hyperbolic horn response character- 
istics are given in Refs. 1 and 2. These expressions are 
derived for frequencies above the cutoff frequency; how- 
ever, these expressions must be modified to correctly 
describe the characteristics at the cutoff frequency and 
below the cutoff frequency. These response character- 
istics are given in the form of impedance expressions; 
therefore, the effect of the horn on the transmission of 
acoustic energy can be described in terms of a lumped 
impedance in an equivalent electrical analog circuit. 

A. Summary of Equations for z T 

Derivations of the required expressions are included in 
Appendix A. Starting with the above cutoff expression for 
the acoustic impedance z T of the general hyperbolic horn, 
the following expressions for z T are obtained: 

(1) For the general hyperbolic horn. 

(a) z T at cutoff. 

(b) z T below cutoff. 

(2) For the special case T = 1, hyperbolic horn (the 

exponential horn). 

(a) z T above cutoff. 

(b) z T at cutoff. 

(c) z r below cutoff. 

The following is a summary of these expressions for z T , 
including the above cutoff expression for the general 
hyperbolic horn. 

1. General hyperbolic horn 
Above cutoff (/ > f c ) : 


\y — tanh (NL) J sin ( bL ) + z M [1 + T tanh (NL)] cos ( bL ) 


+ (— t “ h sin (hL> + [f + { 1 


t 

N(l T l^L| tanh ( NL )J cos ( foL )| 


jkT 


( 1 ) 


At cutoff (/ = f c ): 


Zt 


= ( N [— LAfej, {T + tanh(NL)} + {1 + rtanh(ML)}] +/[NL{ 1 + T tanh (NL) } ] \ 

\ [LN*T{1 + T tanh (NL)} + N {1 + Ttanh (NL)}] + j [LN 2 T 2 z M + {LN 2 Tz u + N (1 - T 2 )} tanh (NL)] / 


( 2 ) 


2 
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Below cutoff (f < f c ): 


f [fc{l + Ttanh (NL)} sinh ( BL )] 

( [{(— B 2 + N 2 T 2 )z„ + k 2 Tz u tanh (NL)} sinh(BL) + {BN(1 - P) Zjr tanh(NL)} cosh(BL)] 

+ j [Nz„ {T + tanh (NL)} sinh (BL) -(Bz M ) {1 + Ttanh (NL)} cosh (BL)}) 

—j [NkT {1 + T tanh (NL)} sinh (BL) + Bfc {1 + Ttanh(NL)} cosh(BL)] J 


2. Special hyperbolic horn (T = 1, exponential horn) 
Above cutoff: 


Zr 


At cutoff: 


_ ( Zu [cos (bL + 0)] + / [ sin (bL)] ~) 
\jz u [sin (bL)] + [cos ( bL — 0)] / 


( 4 ) 


Zr = 




( 5 ) 


horn analysis of paragraph II and the resulting expres- 
sions for the acoustic impedance z T at the throat of a 
given horn. 

The physical parameters of a particular hyperbolic horn 
are determined from the following equation: 

S x = S T [cosh (Nx) + T sinh (Nx)] 2 (8) 

Equation (8) may be written in the exponential form as 
"(e» x + e~* x ) (e* x - e-**)~| 2 


= S T 


+ T- 
2 2 


( 9 ) 


Below cutoff: 


Zr — 


where 


cosh (Nx) = 


e Nx + e -Nz 


sinh (BL) - j i B cosh (BL) - N sinh (BL)} j 

sinh (BL) - _! B »^ ( B c° sh ( BL ) + N sinh ( BL ))]j 

( 6 ) 


All the expressions for z T can be written in the form: 


Zr (<o) — Rr (<o) 4" j Xr (<i>) 


( 7 ) 


where 


Rt (o>) = resistive component of impedance, and repre- 
sents energy transfer by horn 

and 

X r («) — reactive component of acoustic impedance. 

III. Horn Design 

Horns with particular response characteristics can now 
be designed, using the results of the previous hyperbolic 


r>NX g-Nx 

sinh (Nx) = ^ 


Equation (9) may be rewritten 

e y *(l + T) + e- N * (1 - T) - 2 (S*/S T )^ = 0 (10) 

or 

e® (1 + T) ~ 2 (St/Sr)* e v * + (1 - T) = P (x) = 0 

( 11 ) 

the length L of a particular hyperbolic horn is determined 
from Eq. (11) as follows: 

P(*)..l = e iNL (l + T) -2 (S M /S T )*e™ + (1 - T) = 0 

( 12 ) 

Note: S L = Sj^Mouth) 
letting 

e NL = \J/ 
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Therefore, 

L = (1/N) loge ^ 

where ^ is found from Eq. (10), written as: 

r (1 + T) — 2 (S„/S T )* t + (1 - T) = 0 
Therefore, 


* = 


(S i ,/Sr)^±[(S i ,/Sr)-(l-r»)] t 
(1 + T) 


(13) 


(14) 


(15) 


For the physical case of the length L of the horn, the 
root is used: 

, (Su/Sr)* + [Sm/St) - (1 - T*)]V2 

* (ITT) (16) 

When the cross sectional area at the throat S T of the horn 
and the cross sectional area at the mouth of the horn are 
defined, the length L is then only a function of T and N, or 


L = 


where 


r ^ 1 (D + [D 2 - (1 - T*)]V\ 

L ( N , T) = (1/N) loge | (ITT) } 


(17) 


(3) The square root of ratio of the cross sectional area at 
the mouth of the horn to the cross sectional area 
at the throat (S M /S r )\ 

(4) The shape parameter T. 

Figure 1 shows an acoustic horn designed to these speci- 
fications. The equations, shown in Fig. 1, used in the 
design of the horn are: 

y (. x ) — y T [cosh (Nx) 4 - T sinh ( Nx )] (20) 

y (x) — — t/ T [cosh (Nx) + T sinh (Nx)] (21) 



and are derived from the following: 

S x - Sr [cosh (Nx) + T sinh (Nx)] 2 (22) 


D = ( S M /Sr )% 

The value of N is determined from the cutoff frequency 
f c of the horn, using the following expression: 

N = ^ ( 18 ) 

Equation (17) may now be written, in terms of the cutoff 
frequency, as: 

r / c \_ (D + [D 2 — (1 — T 2 )]^) 

L \2nf c ) ° g '' { (1 + T) } (19) 

Equation (19) is, therefore, the expression for the length 
L of a given horn in terms of: 

(1) The cutoff frequency f c% 

(2) The speed of sound c . 


where 

! 7 rt/ 2 (x) (for circular cross section) 

4 y 2 (x) (for square cross section) 
yy 2 (*) (for generalized cross section) 

! 7 rt/| (for circular cross section) 

4t/£ (for square cross section) 

(for generalized cross section) 

! ti ryi (for circular cross section) 

4 yi (for square cross section) 
y yl (for generalized cross section) 

Therefore, 

y y 2 (x) — yt/r [cosh (Nx) + T sinh (Nx)] 2 (23) 

or 

y (x) = ±t/r [cosh (Nx) + T sinh (Nx)] (24) 


4 
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The hyperbolic horn data used for this parametric 
study are given in Table 1. The digital computer program 
written to perform this study using this data is given in 
Appendix B. This program is written for one particular 
group of horns in which the cutoff frequencies f c were 
25, 50, and 100 Hz. Typical curves of the predicted horn 
responses are shown in Fig. 2. The produced curves 
graphically illustrate the computed response of the resis- 
tive component R (o>) for the acoustic impedance versus 
frequency f in relation to the hyperbolic horn family 
parameter T and a fixed ratio of (S M /S r y^. This fixed 
ratio is: 


/ S* V* _ / 49.00 in. 2 Yt 
v Sr ) \1.563 in. 2 / 


D= 5.585. 


(25) 


Table 1. Hyperbolic horn data used in parametric study 


Horn 

T 

L (T) — m 

Horn 

T 

LlD-m 

1 

0.001 

1.3162 

14 

3.00 

0.5944 

2 

0.010 

1.3113 

15 

4.00 

0.4964 

3 

0.100 

1.2646 

16 

5.00 

0.4239 

4 

0.200 

1.2171 

17 

6.00 

0.3685 

5 

0.300 

1.1735 

18 

7.00 

0.3250 

6 

0.400 

1.1333 

19 

8.00 

0.2902 

7 

0.500 

1.0959 

20 

9.00 

0.2617 

8 

0.600 

1.0611 

21 

10.00 

0.2382 

9 

0.700 

1.0285 

22 

20.00 

0.1238 

10 

0.800 

0.9978 

23 

50.00 

0.0501 

11 

0.900 

0.9690 

24 

100.00 

0.0251 

12 

1.00 

0.9417 

25 

500.00 

0.0050 

13 

2.00 

0.7325 

26 

1000.00 

0.0025 




27 

10,000.00 

0.0003 


5.0 

4.5 

4.0 
3 .5 

3.0 

2.5 

1 2.0 

I 

1.5 
1.0 
0.5 

5 0 

(fc* 5.0 

4.5 

4.0 

3.5 

3.0 

2.5 

2.0 
1.5 
1.0 

0.5 
0 

Fig. 2. f-parameter study of hyperbolic horns (frequency response characteristics) 
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The cutoff frequency f c is held constant at 100 Hz, Equation (23) may be written as a function of only L as: 
therefore, 

L = HT) = 0.547, 

(27) 

A P c termination Z M at the mouth of the horn is used 
for the study. This termination is necessary to reveal the 
characteristics of the horn, independent of its particular 
termination. Naturally, the nature of Z M (i.e., zero baffle, 
infinite baffle, etc.) greatly affects Z T for a particular 
configuration. 



[■ P 2 - (1 " 
(1 + T) 



200 400 600 800 1000 1200 1400 1600 1600 2000 'o 200 400 600 800 1000 1200 1400 1600 1800 2000 


FREQUENCY, Hz 

Fig. 2 (contd) 
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Fig. 2 (contd) 
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Fig. 2 (contd) 

A study of the graphic responses in Fig. 2 reveals a with the additional advantage from the averaging effect 

resonant response characteristic that is just above the cut- of several similar, but not identical, horns. This averaging 

off frequency f c which varies with T. With this informa- technique consists of shifting the cutoff frequencies of the 

tion available, acoustic horns can be designed to use this various horns so that the response peaks of some horns 

resonance characteristic. For example, horns with large occur at the response nulls of the remaining horns. This 

resonance characteristics (T < < 1) can be used to assist in effect tends to neutralize the extreme peak-to-null fre- 

simulating acoustic energy spectrums produced by the en- quency spectrum oscillations, 

gines of spacecraft launch vehicles. Non-dimensionalized 

curves of rocket engine noise spectrums are contained in The individual and average frequency response char- 
Ref . 3 and in other discussions on rocket noise. acteristics of four multiple horn arrays, using the horn 

design data listed in Table 2, are graphically portrayed 
in Figs. 3-6, respectively. An overlay example of the 
. esign of Multiple Horn Arrays response characteristics for one array is provided as part 

A multiple horn array can be designed with the desired of Fig. 3. The heavy line in the figure shows the averag- 

resonance response characteristics of a single horn and ing effect resulting from the different horn responses with 
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1.6 


1.4 

1.2 

1.0 

0.8 

0.6 

0.4 

0.2 


3 0 



FREQUENCY, Hz 


Fig. 3. individual and average frequency response characteristics of the horns of array 1 


peaks (resonances) and nulls (relative antiresonances) at 
different frequencies. 


V. Results 

This acoustic horn analysis indicates that the response 
characteristics of a horn can be greatly improved in the 
frequency range just above cutoff by carefully controlling 
the design parameters. This resonance response can vary 
over a large range and consists of several narrow spikes 
which would result in a very rapidly changing narrow 
band analysis. This narrow spike condition can be mini- 
mized by the use of multiple horn arrays. Each horn of 
these arrays must be carefully designed so that the trans- 


fer function of the array will provide the optimum re- 
sponse characteristics. 

VI. Conclusions 

The analysis performed in this study indicates that the 
concept of designing single horns and multiple horn 
arrays to achieve specialized response characteristics 
shows real promise. The analysis assumes linear acoustic 
response, however, this assumption becomes less valid as 
the acoustic levels become very large (above 135 dB in 
reference to 0.0002 jubar). These non-linear effects result 
in harmonic generation which, in turn, affects the spec- 
trums at the higher frequencies. However, these non- 
linear effects have little affect on the resonance peak of 


12 
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FREQUENCY, Hz 


Fig. 3 (contd) 

the acoustic spectrum that is transmitted as a result of 
the horn response characteristics. Also, the horn analysis 
does not allow for the damping effects of either the horn 
or the transmitting medium (air). The effect of such damp- 
ing would decrease the peak amplitudes of the horn 
response function, and therefore, would decrease the 
response functions peak-to-null ratio. Introduction of 
damping could also tend to put a small frequency shift 
in the response function. 

The choice of the parameters for each of the horns of 
the multiple horn arrays is somewhat arbitrary, and is 


Table 2. Design data for multiple horn arrays 


Array 

Horn 

T 

fc , Hz 

L , m 

1 

1 

1.0 

133.3 

0.7065 


2 

1.0 

100.0 

0.9417 


3 

1.0 

72.5 

1.2989 

2 

1 

1.0 

133.3 

0.7065 


2 

1.0 

100.0 

0.9417 


3 

0.5 

72.5 

1.5116 

3 

1 

1.0 

100.0 

0.9417 


2 

1.0 

117.0 

0.8049 


3 

1.0 

134.0 

0.7028 


4 

1.0 

151.0 

0.6237 

4 

1 

1.0 

100.0 

0.9417 


2 

1.0 

114.0 

0.8261 


3 

1.0 

127.0 

0.7415 


4 

1.0 

140.0 

0.6727 


5 

1.0 

153.0 

0.6155 


mainly intended to indicate the general concept of these 
arrays. Future investigations into this concept will include 
developing techniques to determine the design parameters 
for each of the horns in the array so that their overall 
average response function is optimized. The investigation 
will include obtaining the desired resonance response 
function with the smoothest energy spectrum envelope 
possible. 

Future investigations also include verifying these spec- 
trum shaping techniques using a Low-Frequency Plane- 
Wave Sound Generator and Impedance-Measuring De- 
vice, designed and built for this same supporting research 
and advance development task. A report on this device is 
being written and will be published upon final checkout 
of the unit. 
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Appendix A 

Derivation of Horn Response Equations 

I. Derivations 

In this paper the term hyperbolic horn is used to designate those horns whose area law is given by (Refs. 1, 2, and 4) : 

S* = Sr [cosh (Nx) + T sinh (Nx ) ] 2 (A-l) 

The result for the expression of the impedance at the throat z T of a hyperbolic horn as a function of the impedance at 
the mouth z u is given as follows (Ref. 1): 


Zt 


f [ 

'/ jk - NTz 
\ b 

! ) 

+ 1 

nkr _Nz E y 3Lnh(NL) 

sin (bL) +z M [l + Ttanh(NL)] cos (bL) 


(n-I 

(b 2 + N 2 T 2 \ ) 

V jkT ) Z "f + 

( NI '; fe ),a„h(NL)] ! 

iin (bL) + 

ri 

j + < 

( N(l-T*)z u ) 

1 jkT j 

> tanh (NL) 

cos (bL) 


(A-2) 


This equation is valid above cutoff, that is, for k 2 > N 2 (f > / c ). 

Equation (A-2) is rearranged to arrive at expressions for the at-cutoff and below-cutoff frequencies of the general 
hyperbolic horn through derivations Nos. 1 and 2, respectively. These expressions are further simplified for the above- 
cutoff, at-cutoff, and below-cutoff frequencies of the exponential horn through derivations Nos. 2, 3, and 4, respectively. 
Derivations Nos. 5 and 6 are alternate methods of arriving at the at-cutoff and below-cutoff expressions for the expo- 
nential horn. 


A. Derivation No. 1 

At the cutoff frequency, the expression for z T ^> (0 + /0)/(0 + ;0) is an indeterminate form, therefore, L’Hospital Rule 
must be used in determining the value for this expression, as shown by the following derivation. 


Given Eq. (A-2) at b = 0 (cutoff) and z M -> (0 + ;0)/(0 + ;0), thus, by L’Hospital Rule and with k = N; 


(Zt )& = o — 

1 ( ^f" ((— — |~~^) + (-- -fr - ^ Z *) tanh (NL)| sin (bL) + z„ {1 + T tanh (NL)} cos ( bL ) "| ) 

6 = 0 

(A-3) 
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(z T )b=o — j, 


( L [(jk — NTzk) + (/'fcT — Nz^) tanh (NL)] cos (bL) 

\ L 

r (b 2 + N 2 T 2 ) 

L l W ) 

> Zm 

+ [(NT + jkz M ) tanh(NL)] cos (bL) + sin(foL) 

1! 

1 

to 

O - 

N 

* 

L i kT J 


+ [z M (1 + T tanh (NL)} cos (bL) — bLz M {sin 

(bL)} - 

- (bLz M ) {T tanh (NL) sin (bL)}] 


+ 

n r n (i — t 2 ) z M ) 

It \ jkT j 

[ tanh (NL) 

cos (bL) — L 

|>1 

( N (1 — T 2 ) Zj,) 

1 ikT J 

> tanh (NL) 

sin (bL) 


(Zt ) b — 0 rji 


( 

L [(jk - 

- NTz m ) + (jkT - Nz u ) tanh (NL)] +z u [l + T tanh (NL)] 

) r 

K 

(» Mr: 

| + (NT + jkz M ) tanh (NL) 

+ 

J + 

f, N(1 - T 2 )z„) 

r J 

> tanh (NL) 


L [( jk - NTz m ) + (jkT - N Zjf )] tanh (NL) + z M [ 1 + rtanh(NL)] 


]! 

\ 


(zt)(, = „ k | L [iV ( fcT + ?T2NZji/ ) + fcT(NT + jkz M ) tanh (NL)] + [k + (kT + /N( 1 - T 2 )} tanh (NL)]/ 
[ — LNzjf {T + tanh (NL)} + z„ {1 + T tanh (NL)}] + ;[fcL{l + T tanh (NL)}] 


(Z T ) 


/i 


(A-4a) 

(A-4b) 


(A-4c) 

>} 


T)b=0 *\[ LNkT {i + r tanh (NL)} +k{ 1 + Ttanh(NL)}] + / [LN a r-z„ + (Lk 2 Tz„ + N(1 - T 2 )} tanh (NL)}) 

(A-4d) 


Therefore, at f — f c : 

_ ( [-LNz u {T + tanh(NL)} + z M {l + Ttanh(NL)}] +/[NL{1 + Ttanh(NL)}] ) 

Zt “ N \[LN 2 T{ 1 + Ttanh (NL)} + N {1 + Ttanh (NL)}] + / [LN*T*z* + {LN 2 Tz u + N(1 - T 2 )} tanh (NL)] / 

(A-5) 


B. Derivation No. 2 

Below the cutoff frequency, the following derivation is necessary, since for fc 2 < N 2 , b = j ( | b \ 2 )V\ and the various 
terms, such as cos ( bL ), need to be transformed to usable functions. 

Given: 

z T = kX 

( [{(jk - NTzm) + (jkT - Nz m ) tanh (NL)} sin (bL) + bz u {1 + T tanh (NL)} cos (bL)] ) 

\[{NkT + j(b + NT) z M } + {( NkT 2 + jkTz M ) tanh (NL)} sin (bL) ■+ b {k + [kT + jN (1 - T 2 ) Zjf ] } tanh (NL) cos (bL)]f 

(A-6) 


NOTE: Eq. (A-6) is a slightly rearranged version of Eq. (A-2). 
and 


b = jB for fc < N jB = (fc 2 — N 2 ) 1 ^ or B 2 = N 2 — fc 2 sin (jB) = j sinh B cos (jB) = cosh B 


thus 


£r 


= lt 


/ {(jk - NTzm) ± (JkT - Nz^) tanh (NL)} /sinh (BL) 

1 {NfcT + j(-B 2 + N 2 T 2 )z*} + {(NfcT 2 + ;'fe 2 Tz Jf ) tanh (NL)} / sinh (BL) 

+ jBz M {1 + T tanh (NL)} cosh (BL) / 

+ jB {k + + ;'N (1 - T 2 ) z„] tanh (NL)} cosh (BL) / 


(A-7) 
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or 


[k{l + T tanh ( NL )} sinh ( BL )] 

[{(N 2 T 2 - B 2 )z m + k 2 Tz u tanh (NL)} sinh (BL) + {BN ( 1 - T 2 )z M tanh (NL) cosh (BL)}] 

+ i [Nzu {T + tanh ( NL )} sinh (BL) — Bz u {1 + T tanh (NL)} cosh (BL)] | 

- / [NkT (1 + T tanh (NL)} sinh (BL) + Bk{ 1 + Ttanh (NL)} cosh (BL)]/ 


Therefore, for f <f c : 

{ [fc{l + T tanh (NL)} sinh (BL)] 

Zr ~ k \ [{(-B 2 + N 2 T 2 ) z M + k 2 Tz M tanh (NL)} sinh (BL) + {BN(1 - T 2 ) z* tanh (NL)} cosh (BL)] 

+ j [(Nz u ) {T + tanh (NL)} sinh (BL)- (Bz u ) {1 + Ttanh (NL)} cosh (BL)]} 

- / [NfcT {1 + Ttanh(NL)} sinh(BL) + Bk {1 + Ttanh (NL)} cosh(BL)]/ 


The expression of the impedance at the throat z T for an exponential horn (T = 1), in terms of the impedance at the 
mouth z u , is given (Ref. 2) as follows (above cutoff) : 


/ pc\ ( SgZjf [cos (BL + 9)] + / (pc) (sin (BL)]) 
\ Sj / ) jS 2 Z„ [sin ( bL )] + ( P c) [cos (bL — 0)] / 


assuming 


Therefore, 


Zt = 


Z T Zji 

Zt ~7^7sJ Zm ~T^sT) 


z u [cos (bL + 9)] + j [sin (BL)]/ 
jz u [sin (bL)] + [cos (bL — 0)] / 


(A-ll) 


C. Derivation No. 3 

Since at the cutoff frequency z T -> (0 + /0)/(0 + jO), L’Hospital Rule is applied for the following result. 
Given Eq. (A-ll) at the cutoff frequency f c (from 4k 2 = M 2 ) with 


b = y 2 (4k 2 - M 2 )K = (k 2 - N 2 )^ = 0 0 = tarn 5 (N/b) = tt/2 

Therefore, at f = f c : 

_ ( Zu [cos (tt/2)] + / [sin (0)] ) 0 + ;0 

Zt \jz u [sin(O)] + [cos (-tt/2)]} 0 + f0 


(A-12) 
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Hence, by L’Hospital Rule, z T (at f = f c ) is found as follows: 

r (d/dbKM [cos (bL + fl)] + / [sin (bL)]} ~l 
Zt |_ (d/db)(jz u [sin (bL)] + [cos (bL — «)]} J 6=0 

_ f (d/db){zu [cos (bL + tan' 1 N/b)] + j [sin (bL)] } ~l 
“ |_(d/dB){/z* [sin (bL)] + [cos (bL - tan" 1 N/b)] } J b=0 

i-z u [sin (bL + tan- 1 N/b)] j^L + jj- — J + / ( L ) [cos ( bL )]^ 

| f(x K L) [cos (bL)] - [sin (bL -tan- 1 N/b)] [ L + (T+^NyP)] 

j-z u [sin (bL + tan" 1 N/b)] £ - + ^ J + /(L) [cos (bL)]^ 

— \ r l (tf 2 _j_ £ 2 ) _|_ ^“1 / 

^ /(z*L) [cos (foL)] - [sin (bL — tan -1 N/b)] + ^ J ) ^ # 

~(z u ) [sin (tt/ 2)] ( LN n 7 N ) + /(*•) [cos (0)] 

/(zjtfL) [cos (0)] - sin (-tt/2) 

-M(L-±)+m +fa)(^-L) + i(L) ^(l- J f 1 ) + <( j T 1 ) 

j(,„L) + (l + i) HzM + (l + -^) i*»( i f-) + ( tl t + 1 ) 

Therefore, at f = f c : 

/ LM\ , ,/ML\ 

2 ) + , V 2 ) 

Zt ~ / ML\ . / LAf\ 

( 1 + — ) + >*'(-r) 

D. Derivation No. 4 

Below the cutoff frequency, that is for 2fc < M(or fc < N), following results are derived: 

Given Eq. (A-ll) and the following definitions of trigonometric identities: 

(a) b = jB (b) cos (8 ± jB) = cos (6) cosh (B) q: j sin (0) sinh (B) 

(c) sin (;'B) = / sinh B (d) cos (jB) = cosh B 

(e) cos (±8 + bL) = cos (±0 + jBL) = cos (8) cosh (BL) qi j sin (8) sinh (B) 

(f) 8 = tan- 1 (N/b) = tarn 1 (N/jB) (g> cos (8) = jB/(N* - B 2 ) l/2 (h) sin (0) = N/(N 2 - B 2 ) 1 ^ 


(A-13a) 

(A-13b) 

(A-13c) 

(A-13d) 

(A-13e) 

(A-13f) 

(A-14) 
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then 


__ f z M [cos (0) cosh (BL) — / sin (0) sinh ( BL )] + y [y sinh ( BL )] ^ 
Zt \jZM [/sinh(BL)] + [cos (0) cosh (BL) + / (sin (0) sinh (BL)}]J 

Therefore, for / < fc- 


( [sinh (BL) - /T ft, _* {B cosh (BL) - N sinh (BL)}] ) 

Zr = < 1 P > (A-16) 

([2* sinh (BL) - j - _ (B cosh (BL) + N sinh (BL)}]) 


E. Derivation No. 5 

Equation (A-14) can also be derived from the general expression for z r using the family parameter T equal to 1. Given 
Eq. (A-5) with T = 1 for the exponential horn, then Eq. (A-5) reduces to: 


(1-LNzm (1 + tanh (NL)} + z„ { 1 + tanh(NL)}] + / [NL {1 + tanh (NL)}]) 
ZT ~ k \ [ LN 2 {1 + tanh (NL)} + N{ 1 + tanh (NL)}] + j [LN 2 z„ {1 + tanh(NL)}] ) 

z u ( 1 - LN) + / (NL) 

Zt (1 + LN) + jz„(LN) 

however, N = M/2; therefore. 


(A-17) 


(A-18) 


Zr — 


/ LM\ 

i ^ f ML \ 

z “ y 1 2~ ) 

i +, hr) 




Eq. (A-14) 


F. Derivation No. 6 

Equation (A-16) can also be derived from the general expression for z T using the family parameter T equal to 1. Given 
Eq. (A-9), then, for T = 1: 


( [fc sinh (BL)] + / [z u { N sinh (BL) - B cosh (BL)}] ( 
\[ Zjf fc 2 sinh (BL)] - / [fc (N sinh (BL) + B cosh (BL)}] j 


for /B — b or B = — jb : 


( [fcsinh(BL)] + / [z^ {N sinh (BL) + /b cosh (BL)}] ) 
Zr _ * t [Ziffc 2 sinh (BL)] - f [fc {N sinh (BL) - jb cosh (BL)}] J 


(A-20) 


with tan- 1 (N/b) = 0, then 
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( [( k/b ) sinh (BL)] + j [z u ( N/b ) sinh (BL) + / cosh (BL)}] 1 

Zt ~ k \[zMk {(k/b) sinh (BL)}] - j [k {(N/b) sinh (BL) - /cosh(BL)}]/ 

with sec 0 — k/b , then 

— v / [sec (0) sinh (BL)] + ; [a^ {tan (0) sinh (BL) + /cosh (BL)}] ) 

%r {sec (0) sinh (BL)}] — / [fc {tan (0) sinh (BL) — / cosh (BL)}]J 

__ J [sinh (BL)] + / [z^ {sin (0) sinh (BL) + j cos (0) cosh (BL)}] | 

2t \[zjtA: {sinh (BL)}] — jk [sin (0) sinh (BL) — j {cos (0) cosh (BL)}] j 

however, cos 0 = jB/(N 2 — B 2 ) 1 ^ and sin 0 = N/(N 2 — B 2 ) 1 ^; therefore, 


(A-21) 


(A-22) 

(A-23) 


( [sinh (BL) - / _* {Bcosh(BL) - N sinh (BL)}] j 

( > - Eq. (A-16) 

([z M smh(BL) - {Bcosh(BL) + N sinh (BL)}]) 
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Appendix B 

Digital Computer Program 


c 

X- 



ER000300 

c 


ACOUSTICAL HORN PROBLEM APRIL, 1966 

xxxxxxxx 

COMPLEX T.N,L,F»B,Tl»IC»A»ZM,T2»T3.SINH»COSH,TANH»T4»T5»T6»T7. 
1ZT,ZT1,C,SMALLB, SINE, COSINE, T8, T9, T10, X. X2. THETA, NUM 

c 

c 


THESE ARE CONSTANTS USED THROUGHOUT THE PROGRAM 

ER000200 

C= ( 344.0,0.0 
A= (0. 10,0.0 ) 

c 

c 



ER000500 


101 

IFIRBT=0 


C 


READ IN VARIABLES - HORN NUMBER, T, N AND L. 
READ (5,1) I HORN , T, N, L 


C 

1 

FORMAT (15. 6 F 1 0 * 4 ) 

— -** 


c 


IF HORN NUMBER IS EQUAL TO ZERO, NO MORE DATA. 
IF ( I HORN ) 1000,1000,2 


c 


PRINT OUT HEADER 


C 


WRITF (6,4) IHORN,T,N,L 


c 


V 

FORMAT MH1, IPX. 21 HACOUST I C HORN PROBLEM, / , 2 8X , 4HH0RN , 


_c_ 


112, 5H T = , 2F5.2* 5H N = , 2F5.2, 5H L = , 2F5.2) 

xxxxxxxx 

c 


IF HORN IS ONE, TWO, OR THREE, DO FORMULAS 9, 3, AND 2 
GO TO (6,6,6,100,100) ,IHORN 


c 

c 

c 

— - A 

SET F TO 5 CPS, GO TO EQUATION 9 
F=(S. 0.0.0) 



402 

T 1 =N*L 

T2 = CFXP(T1 1 


T 3 = CEXP ( -T 1 > 

TANH= ( T7-T31 J ( T? + T3J_ 

C 

_c_ 



I NO 1 2 8 00 

c 

? 

IF F HAS REACHED 100 CPS, GO TO EQUATION 3 

FF = RFAL(F) 


IL. 


IF (FF-100.0) 8,20,20 

IN012600 

C 


COMPUTE B, ZM, TANH, COSH, AND SINH 



r 

8 B= ( 1. 83,0. 0)*C SORT ( (1.0,0.0)-F*F*(0. 0001, 0.0) ) 


c 

COMPUTE K 

. v - i . m «•» ,n . a > #f 


JO 

T1=K*A 

INQ-1 1600 

ZM=( 0.24,0. 0)*Tl*Tl + < 0.0, 1.0)*<0.57,0.0)*T1 

Ti-n*i 

T 2=CEXP ( T 1 ) 

t a = rr yd / -T 1 \ 

SINH=(T2-T3)/(2. 0,0.0) 

CO SH=(T2+T3)/(2. 0,0.0) 
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c IN011400 

c 


c 

COMPUTF FOUATION 9 



o T1=K*( ( (1.0,0. 0)+T*TANH)*SINH) 
T? = ( N*ZM* ( T+TANH ) *S I NH ) 

T3=-R*ZM*( ( I .0,0.0)+T*TANH)*COSH 
T 4=( (-B*R+N*N*T*T I #ZM+K*K*T *Z M*T ANH ) *S I N H 

T 5 = { B*N* ( (1.0,0. )-T*T)*ZM*TANH)*COSH 

T A=N*K*T* ( M.0,0.0)+T»TANH)»STNH _ 

f 

T7 = B#K* ( ( 1.0*0«0)+T*TANH)*COSH 


I NOO 9 5 00 

c 



IN007215 

C 

FORM 7 T 


I N009200 


ZT=K* ( (T1+(0.0,1.0)*(T2+T3) )/(T4+T5-(0.0»l 

,0)*( T6+T7) ) ) 

IN007205 

C 

IS THIS THE FIRST TIME THROUGH 
IF f IF I RST ) 14, 14., 16 



c 

YES - GO BACK, COMPUTE EQUATION WITH OTHER 
14 TFIRST = 1 

ZM 


ZT 1 =ZT 

ZM= ( 1.0, 0.0) 

JC_ 

GO TO 9 


I N007500 


C I N007300 

C PRINT OUT BOTH VA LUES WITH F 

1A FF=REAL(F) ' 

WRITE (6 .18) FF.ZT1.ZT 

18 FORMAT (1H0, 2X , F6.1.4H CPS, 3X, 2E17.5, 3X.2E17.5) 

C ADD 5 CPS TO F 

F=F+(5. 0,0.0) 

I F I RST = 0 

GO TO 7 


C 

£ 

C #- 

£ COMPUTE FQUATION 3 


C FIRST RESET ZM 

£ COMPUTE X 



20 

<=(.0183, 0.0 ) *F 



_2_L 

__T 1 = K*A 




ZM=(0.24, 0.0)*T1*T 1+10.0, 1.0) *(0.5 7, 0. 0 )*T1 

IN003500 

c 

__Z2_ 

T 1 3 -L*N*ZM# ( T+TANH) 

I N003300 

T2=ZM*( (1.0,0.0)+T*TANH) 
T3=N*L*< ( 1.0,0.0)+T*TANH) 

T4=L*N*N*T*( (1.0,0.0)+T*TANH) 
T5=N*( ( 1.0,0. 0)+T*TANH) 

T6=L*N*N*T*T*ZM 

T7=(L*N*N*T*ZM+N« ( ( 1.0,0.0)-T*T ) ) *T ANH 


ZT = N*( (T1+T2+(0.0,1.0)*T3 > / ( T4 + T 5+ ( 0 . 0 , 1 . 0 ) * ( T6 + T 7 ) ) ) 


IN007100 
I N006900 


24 
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xxxxxxxx 


c 



XXXXXXXX 

IN000800 

c 


SUBSTITUTE OTHER VALUE OF ZM 
LF (IFIRST! 24,24j2_6 . 



24 

IFIRST=1 

ZT 1 = 7 T .... _ _ .... ..... 


2M=(1. 0,0.0) 

GO TO 7 7 . . . . . 

c 

_26 

FF = RFA! (F) ..... _ . ... . 

IN000200 

ifir$t=o 

WRITF nS.18)FF.ZT1.2T _ 

c 

c_ 


NOW COMPUTE FQUAT ION 2 

XXXXXXXX 

c 

28 

DETERMINE IF F IS EQUAL TO 150 CPS 
EF=REAL(F) 




IF (FF-150.0) 30,32,32 

I N000300 

c 

80 

F UNDER 150 CPS - ADD 2 CPS 
F=F+( 2.C.0.0) 


£_ 


GO TO 36 


C 

8? 

ADD 5 CPS TO F 
E=F+(5. 0,0.0) 


C 



XXXXXXXX 


C 

.. 34 

DETERMINE IF F = 500 CPS 

FF=R£4L(P) 


jC_ 


IF (FF-505.0) 36,50,50 


c 

.C— 


rnMDMTtr rnMATTAW 7 f,PT F T P F T \/AflfF OF 7 M 



36 

K= (. 0183*0. 0)*F 

T 1 =k-*4 



— 3.7 

ZM= ( 0.24,0.0 ) *Tl*Tl+( 0.0, 1.0) »( 0.57 *0.0 )*T1 

SMAI l R- I 1 .ftS.0.0 l*rSQRTI { 0.0001 .0.0)*F*F-( 1 .0,0.0 ) ) 


SINE=CSIN(SMALLB*L> 
cnsjNF=rrosr smai l r*l ) 


38 

T 1 = ( -N*T#ZM+ ( 0 • 0 » 1 • 0 ) *K ) /SMALLB 
-T?=f-N*7M+( 0.0.1 .0 )*T*< )/ SMALL B ... 


T3=(T1+T2*TANH)*SINE 

T4=7M#I M .O.fi.O)+T*T ANH ) *005 T NF 

T5 = ( SMALLB*SMALLB+N*N*T *T ) / ( ( .0 , 1 .0 ) *K*T > 
TA=fN-T5*7M) /SMAII R 

T7=(N*T+( 0.0,1 .0)*K*ZM)/SMALLB 
Ta= IT6+T7*TANH)*STNF 

T9=(1.C,0.0)-(N*(( 1.0, 0.0 )-T*T)*ZM)/ ( <0.0,1.0)*K*T) 
.1X0 = f ( 1 .0,0.0) /T+T9*TANH)*C0SINE 

_C_ 


ZT=(T3+T4)/(T*(T8+T10)) 


c 

x_ 


GO rnMPUTF FOIIATTON AGAIN FOR SFCOND ZM 



-47 

IF (IFIRST) 42*42,44 

I F I RST = 1 ... . . .. _ 


2M=(1. 0,0.0) 

_ ZT 1 =7 T . 

GO TO 38 
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WRITE (6.18IFF.ZT1 ,ZT 

tftpst=Q — 

GO TO 28 

-»» 

IF THIS IS HORN 3. GO THROUGH 2ND SET OF EQUATIONS 

50 GO TO (101*101 *51 ) .IHORN 

PRIMT SEPARATOR 

51 WRITE (6. 52) 

5? FQRMAT(50H ) 


FIRST SET F TO 5 CPS 


iPUTE ZM 
C. 0183*0. 0)*F 


CC=REAL( (B+N)/(B-N) I 

CC =AI QG 1 CC*CC 1 / 4.0 __ 

T1=CMPI_X(CC,0.0) 

T 8 =CF XP ( T 1 ) 

T9 = CEXP(-T1 ) 

FORM COSH AND SINH OF THE IN 
T2=(T8+T9 ) / ( 2.0 .0.0) 


54 ZT=(ZM*(T2*T6-T3*T7)-T7)/ ( -ZM*T 7+T2*T6+T 3*T7 ) 
IS THIS THE FIRST TIME WITH THIS CPS 


FIRST TIME - SAVE ZT 


ZM=(1. 0,0.0) 


GO TO 54 


WRITE (6,18) FF.ZT1.ZT 

DETERMINE IF CUT OFF HAS BEEN REACHED. 

E = F + ( 5 « 0 .0.01 

I F I RST = 0 

FF=REAL(F) 

IF (FF-100.0) 53.66.66 

COMPUTE EQUATION 17 

K=t. 0183.0, 0)»F 

T1=K*A 

7M= (0.24. 0.0 )»Tl«Tl+( 0.0, 1.0 )»( 0.57,0.0 )»T1 

I FIRST a o 


XXXXXXXX 
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Nomenclature 


a radius of cross sectional area at 
mouth, m 

b + y 2 (4k 2 - Af 2 )^ = + (k 2 - 2V 2 )^, nr 1 

B + (2V 2 — kr) l & = ;fe, m -1 

c velocity of sound in medium trans- 
mitting sound waves (air), m/s 

f frequency of sound waves, Hz 

f c cutoff frequency of sound waves, Hz 

k wave number = 2 tt/A = *>/c, m -1 

L length of horn, m 

M flare constant of exponential horn = 
22V, m -1 

2V flare constant of hyperbolic horn, nr 1 
R t from z T = Rt + jX T , N-s-m -6 


S x cross sectional area of hyperbolic (and 
exponential) horns at a position x units 
from throat, m 2 

S M = cross sectional area at mouth, 
m 2 

S T — cross sectional area at throat, 
m 2 

X T from z T = Rt + /X T , N-s-m -5 

Zi ( i = M, r) acoustical impedances of various 
acoustic systems to be described, Zi = 

%i/(pc/Si) 

z M — acoustical impedance at mouth 
z r = acoustical impedance at throat 


Zi ( i = 2Vf , r) acoustical impedances of various 
acoustic systems to be described (mks 
acoustic Q), N-s-m -5 

= acoustical impedance at mouth 
of particular horn being de- 
scribed, N-s-m -5 

z T = acoustical impedance at throat 
of particular horn being de- 
scribed, N-s-m -5 

£ tan -1 (Xr/Rr), deg 
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T2»(L*(3.66*0.0 ))/ (2.0.0.0) 


IF (IFIRST) 74,74,78 


ZM=( 1 .0,0.0 ) 


GO TO 72 


WRITE (6*18) FFiZTlfZT 


F =F+ ( 2. 0,0.0) 


780 K= ( .0183,0. 0)*F 


ZM= (0. 24*0.0 )*T1*T1+(0.0»1.0)*< 0.57*0.0 )*T1 


5MALLB=( 1 . 83 , 0 . 0 ) *C SORT ( { . 0 00 1 , 0 . 0 ) *F*F - ( 1.0 , 0 . 0 ) ) 


CC=REAL(N/SMALLB) 


THETA=CMPLX (CC ,0.0) 


GOT THETA - NOW COMPUTE EQUATION 16 

FIRST GET SINES AND COSINES 

-2&7 SINE * £S IN 1SMALL B «-U 

COSINE = CCOS(SMALLB*L+THETA) 
Tl=CCOS(SMALLR*l_-THETA ) 


T2=ZM*C0SINE+( 0.0* 1.0 )*SINE 
T 3 = ( 0.0, 1 .0 )»ZM»SINE+T1 
ZT=T2/T3 


GO THRU FOR OTHER ZM 
IF (IFIRST) 788,788,789 


788 IFIRST*! 


ZM= ( 1 .0,0.0) 


GO TO 787 
8 9 FF = REAL(F ) 


I F I RST = 0 

WRITE (6,18) FF,ZT1,ZT 


86 F = F+( 5. 0,0.0) 
GO TO 780 


00 CONTINUE 


1000 STOP 
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